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Abstract. A metrized graph is a compact singular 1-manifold endowed with 
a metric. A given metrized graph can be modelled by a family of weighted 
combinatorial graphs. If one chooses a sequence of models from this family such 
that the vertices become uniformly distributed on the metrized graph, then the 
ith largest eigenvalue of the Laplacian matrices of these combinatorial graphs 
converges to the ith largest eigenvalue of the continuous Laplacian operator on 
the metrized graph upon suitable scaling. The eigenvectors of these matrices 
can be viewed as functions on the metrized graph by linear interpolation. These 
interpolated functions form a normal family, any convergent subsequence of 
which limits to an eigenfunction of the continuous Laplacian operator on the 
metrized graph. 
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1. Introduction 

Roughly speaking, a metrized graph T is a compact singular 1-manifold endowed 
with a metric. A given metrized graph can be modelled by a family of weighted 
combinatorial graphs by marking a finite number of points on the metrized graph 
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and declaring them to be vertices. On each of these models we have a Lapla- 
cian matrix, which acts on functions defined on the vertices of the model. On the 
metrized graph there is a measure-valued Laplacian operator. The goal of this 
paper is to prove that the spectra of these two operators are intimately related. 
Indeed, we will show that if we pick a sequence of models for T whose vertices be- 
come equidistributed, then the eigenvalues of the discrete Laplacians on the models 
converge to the eigenvalues of the Laplacian operator on T provided we scale them 
correctly. Moreover, we can show that in a precise sense the eigenvectors of the 
discrete Laplacian converge uniformly to eigenfunctions of the Laplacian operator 
on the metrized graph. 

This type of approximation result for Laplacian eigenvalues dates back at least 
to the papers of Fukaya ([KeFu]) and Fujiwara ([KoFul], [KoFu2]). In [KeFu] the 
measured Hausdorff topology is defined on closed Ricmannian manifolds of a fixed 
dimension subject to certain curvature hypotheses, and it is shown that convergence 
of manifolds in this topology implies convergence of eigenvalues for the associated 
Laplace-Beltrami operators. An analogue of the measured Hausdorff topology for 
the class of finite weighted graphs is defined in [KoFul]; a similar convergence result 
for eigenvalues of the discrete Laplacian operator is obtained. In [KoFu2], Fujiwara 
approximates a closed Riemannian manifold by embedded finite graphs and proves 
that the eigenvalues of the Laplace-Beltrami operator can be bounded in terms of 
the eigenvalues of the associated graph Laplacians. In the present paper we use an 
approach remarkably similar to the one in [KoFul]; this is purely a coincidence as 
the author was unaware of Fujiwara's work until after giving a proof of the main 
theorem. It would be interesting to sec if the methods employed here can improve 
upon [KoFu2] when applied to the approximation of Riemannian manifolds by 
graphs. 

Without giving too many of the details, let us briefly display some of the con- 
tent of the main theorem (Theorem 1) in a special case. See section 2 for precise 
definitions of all of the objects mentioned here. Let T = [0, 1] be the interval of 
length 1. Let Gat be a weighted graph with vertex set Vn = {<7i, • • • , (In}, edge set 
En — {qiq i+1 : i = 1, . . . , N — 1}, and weight N — 1 on each edge. The reciprocal of 
the weight of an edge will be its length. We say that Gn is a model of our metrized 
graph T; see Figure 1. 

For each N we can define the Laplacian matrix Q N associated to the graph Gn- 
It is an N x N matrix which contains the weight and incidence data for the graph. 
Let A^jv be the smallest positive eigenvalue of Qn- The following table gives the 
value of NXi^n (the scaled eigenvalue) for several choices of N. In each case, the 
eigenspace associated to this eigenvalue has dimension 1. 



N 


NXi t N (approximate) 


5 


7.6393 


10 


8.8098 


50 


9.6690 


100 


9.7701 


200 


9.8201 


500 


9.8498 



Now let / : T — > C be a continuous function that is smooth away from a finite 
set of points and that has one-sided derivatives at all of its singular points. The 
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Laplacian of such a function / is defined to be 

A/ = -f"{x)dx - J2 S p> 

p a singular 
point of / 

where dx is the Lebesgue measure on the interval, S p is the point mass at p, and 
<Tp(,f) is the sum of the one-sided derivatives of / at the singular point p. We say 
that a nonzero function / is an eigenfunction for A with respect to the measure dx 
if the following two conditions obtain: 

• There exists an eigenvalue A > such that A/ = Xf(x)dx. 

The eigenvalues of A with respect to the measure dx are n 2 n 2 for n = 1,2,3,..., 
and the eigenspace associated to each eigenvalue has dimension 1. Denoting the 
smallest eigenvalue by Ai(T), we see that Ai(T) = n 2 w 9.8696. The values of 
NXi.n in the above table could reasonably be converging to Ai(r). 

Part (A) of Theorem 1 asserts that the scaled eigenvalues AA^jv do indeed 
converge to Ai(T). A similar statement can be made for the second smallest eigen- 
values of Qm and A, as well as the third, etc. The fact that the dimensions of the 
eigenspaces for \\.n and Ai(r) agree is no coincidence, and a precise statement of 
this phenomenon constitutes part (B) of Theorem 1. If we choose an ^-normalized 
eigenvector Hn of the matrix Qjy for each N, then the sequence {hjy} can be viewed 
as a family of piecewise affine functions on the interval by linear interpolation. Part 
(C) of the main result states that this family is normal, and any subsequential limit 
of it will be an L 2 -normalized eigenfunction for A with respect to the dx measure. 

We make all of this precise in the next section after defining some of the necessary 
notation and conventions. We will follow quite closely the treatment of metrized 
graphs given in [BR] . Applications of metrized graphs to other areas of mathematics 
and the physical sciences can be found in [Ku] and [BR] . For a more conversational 
approach to metrized graphs, see the expository article [BF] . The reader should be 
aware that metrized graphs also appear in the literature under the names quantum 
graphs, metric graphs, and c 2 -networks. 



[0,1] G 



N 



1/4 



Ql <?2 <?3 <?4 <?5 



FIGURE 1. Here we see the metrized graph T and the graph Gjy 
where N = 5. Each edge of Gn has length 1/4 or weight 4. One 
can view Gat as a discrete approximation to the metric space T. 



2. Definitions, notation, and statement of the main theorem 

A metrized graph T is a compact connected metric space such that for each point 
per there exists a radius r p > and a valence n p e N such that the open ball 
in r of radius r p about p is isometric to the star-shaped set ^ re 2 ^ lm / n p : < r < 
r p , 1 < m < n p } C C endowed with the path metric. A vertex set V for T is any 
finite nonempty subset satisfying the following properties: (i) V contains all points 
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per with n p 2; (ii) for each connected component Ui C T \ V, the closure 
= Ui is isometric to a closed interval (not a circle); and (iii) the intersection 
ei Dej consists of at most one point when i ^ j. The set will be called a segment 
o/r with respect to the vertex set V. Given a vertex set V for T, one can associate 
a combinatorial weighted graph G = G(V) called a model for T. Indeed, index the 
vertices of G by the points in V and connect two vertices p,q'm. G with an edge of 
weight 1/L if r has a segment e of length L with endpoints p, q. The hypotheses 
we have placed on a vertex set ensure that G is a finite connected weighted graph 
with no multiple or loop edges. 

Given a vertex set V for T, a segment of length L can be isometrically parametrized 
by a closed interval [0, L\. This parametrization is unique up to a choice of orien- 
tation, and so there is a well-defined notion of Lebesgue measure on the segment 
with total mass L. Defining the Lebesgue measure for each of the finite number 
of segments gives the Lebesgue measure on F, which we denote by dx. Evidently 
it is independent of the choice of vertex set. For the scope of this paper we will 
assume that T is a fixed metrized graph on which the metric has been scaled so 
that J r dx = 1; i.e., T has total length 1. 

Choose a signed measure of total mass 1 on F of the form 

n 

(1) ii = u}(x)dx + ^c j 5 P] (x), 

where w is a real- valued piecewise continuous function in i 1 (r), ci, . . . , c„ are real 
numbers, and S p (x) denotes the unit point mass at p £ T. We may assume that 
X = {Pi; • • • >Pn} is a vertex set for F containing all points where u) is discontinuous. 
This particular vertex set X will be fixed for the rest of the article. The choice of 
measure ix allows some flexibility in applications of the theory (cf. §14 of [BR]). 

For each p £ T, we define the set Vec(p) of formal unit vectors emanating from 
p. This set has n p elements in it, where n p is the valence of T at p. For v <E Vec(p), 
we write p + ev for the point at distance e from p in the direction v, a notion which 
makes sense for e sufficiently small. Given a function / : T — > C, a point p G T, 
and a direction v e Vec(p), the derivative of f at p in the direction v (or simply 
directional derivative), written D$f(p), is given by 

D 9 f{p)= lim^M, 

provided this limit exists. 

The Zhang space, denoted Zh(r), is the set of all continuous functions / : T — > C 
for which there exists a vertex set Xf for T such that the restriction of / to any 
connected component of T \ Xf is C 2 , and for which f"(x) e i 1 (r). Directional 
derivatives of functions in the Zhang space exist for all points of T and all directions. 
The importance of this space of functions will be made clear in just a moment. 

The Laplacian A(/) of a function / £ Zh(r) is the measure defined by 

A(/) = -f"(x)dx £ £ A,/(p) 1 5 p (x). 

per [iTeVec(p) J 

One sees that X^eVcc(p) D$f(p) = for any point p £ T \ Xf so that the outer sum 
above is actually finite. For any f,g £ Zh(r), the Laplacian operator satisfies the 
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identities 



(2) 




•x. 



Letting g be the constant function with value 1 in this last expression shows that 
the measure A(/) has total mass zero. 

Define Zh M (r) to be the subspace of the Zhang space consisting of all functions 
which are orthogonal to the measure Li; i.e., all / £ Zh(r) such that J T fdfi = 0. 
A nonzero function / £ Zh Ai (r) is called an eigenfunction of the Laplacian (with 
respect to the measure /j.) if there exists an eigenvalue A £ C so that 



The integral on the left side of this equation is called the Dirichlet inner product of 
/ and g and is denoted (/, g) Dir ; the integral on the right side is the usual L 2 -inner 
product and will be denoted (f,g) L2 - Thus we can restate the defining relation 
for an eigenfunction as (f,g} uiT = X(f,g) L2 - Setting g — f and applying the 
relations in (2) shows that each eigenvalue of the Laplacian A must be real and 
positive. The eigenvalues constitute a sequence tending to infinity; in particular, 
the dimension of the eigenspace associated to a given eigenvalue is finite. Let 
< Ai(r) < A2(r) < As(r) < • • • denote the eigenvalues of A with respect to the 
measure /i. The dimension of the eigenspace corresponding to A, (T) will be denoted 
di(T). Even though we have suppressed ll in the notation for the eigenvalues, they 
do depend heavily on the choice of measure. 

The Laplacian operator can be defined on a much larger class of continuous func- 
tions than Zh(r), and the eigenfunctions of A in general need not lie in the Zhang 
space. However, due to the "smoothness" of our measure fx, every eigenfunction of 
A is indeed C 2 on the complement of our fixed vertex set X . In fact, if lo(x) (the 
continuous part of fx) is C m on r \ X , then each eigenfunction of the Laplacian is 
Qm+2 on conl pi cmt ; n t of X . All of this is illuminated by Proposition 15.1 of 



Now we turn out attention to the discrete approximation of T by its models. For 
each positive integer N > #X, choose a vertex set X C Vat C T in such a way 
that #Vat = N. We will add another hypothesis to this vertex set momentarily 
Define the measure dxN on T to be the probability measure with a point mass of 
weight l/N at each point in Vn. Choose the vertex sets {Vn} so that the sequence 
of measures {dxff} converges weakly to the Lebesgue measure dx. In doing this, 
we are modelling the metrized graph T by finite weighted graphs Gn = G(Vn) 
whose vertices become equidistributed in T. This sequence of models will be of 
great interest to us. 

For each N, we also choose a discrete signed measure /in supported on Vn with 
total mass 1 and finite total variation such that the sequence {lin} tends weakly 
to ll. For example, we could construct such a sequence in the following way. For 
each p £ Vn , consider the set of points of T which are closer to p than to any other 
vertex: 

Ap = {x £ T : dist(x,p) < dist(x, q) for all vertices q £ Vn \ {p}} , 

where dist(-, •) is the metric on T. Note that this set is Borel measurable (in fact, it 
is open). We define the discretization of the measure \x associated to Vat by setting 




for each g £ Zh M (r). 



[BR]. 
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Hn{p) = fJ-(Ap) for each p £ Vat. The sets are pairwise disjoint, so this yields 
a discrete signed measure on T supported on the vertex set Vn- As X C Vjv, we 
can assert that [In has total mass 1 (because then all point masses of \i are picked 
up by some Ap). Recalling the definition of the measure \i in (1), we see that the 
total variation of \in satisfies 

M(r) < / \u(x)\dx + 

Jt i= i 

One can check that the measures {hn} tend weakly to /i as N tends to infinity. 

We now wish to tie together the two notions of functions on the vertices of a 
model G and functions on the metrized graph T. The most natural way to do this 
is via linear interpolation of the values of a given function on a vertex set. The class 
of continuous piecewise affine functions on V, denoted CPA(r), is defined to be the 
set of all continuous functions / : T — > C for which there exists a vertex set Xf 
with the property that / is affine on any connected component of the complement 
of Xf. More precisely, if U is a connected component of T \ Xf, then its closure 
e = U admits an isometric parametrization s e : [0, L] — > e. To say that / is affine 
on e is to say that there are complex constants A, B depending on the segment e 
so that / o s e (t) = At + B. For each vertex set V of T, we define Funct(y) to be 
the subclass of CPA(r) whose values are determined by their values on V. That is, 
Funct(y) is any continuous function on T which is a linear interpolation of some set 
of values on the vertex set V. There is a natural identification of Funct(Vjv) with 
the complex vector space C^. We also define the £ 2 -inner product of two elements 
in Funct(V/v) to be (/, g) e2 = J T f(x)g(x)dxN ■ It is important to note which value 
of N we are considering when computing this inner product. 

Each of our preferred models Gat is equipped with a combinatorial weighted 
Laplacian matrix (or Kirchhoff matrix), denoted by Qn, which we can view as an 
abstract linear operator on the space of functions Funct(V/v). Label the vertices of 
Gn by gi,... , qNi and assume that edge qiqj = qjqi has weight Wij (recall that Wij 
is the reciprocal of the length of the segment qiq^ in T and wu = for all i). By 
definition (cf. [Mo] or [Bo]), the entries of Qn are given by 

—Wij, if i 7^ j and is adjacent to qj 
0, if qt is not adjacent to qj. 

For / £ Funct(Viv), we write Qat/ or {QNf}(x) for the unique function in Funct(Viv) 
which has value J2j[QN]ijf{qj) at the vertex qi. Evidently the function Qjv/ ex- 
hibits the same information one gets from multiplying the matrix Qn on the right 
by the column vector with jth entry f(qj). For functions in Funct(Viv), the Lapla- 
cian matrix is closely related to the Laplacian operator A via the formula 

n 

(4) A(f) = J2{QNf}(qj)S qj . 

j'=i 

This is an easy consequence of the definitions of all of the objects involved (cf . [BF] , 
Theorem 4). 



(3) [Q N ]ij = { 
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In analogy with the setup for the continuous Laplacian, we define the class of 
functions Funct^ w (V/v) to be the subclass of Punct(Vjv) orthogonal to the mea- 
sure fijy (a subspace of complex dimension N — 1). That is, / e Funct ([tJV (Vjv) if 
Jp fdfiN = 0. A nonzero function / € Funct MN (V/v) will be called an eigenfunction 
for the discrete Laplacian Qn (with respect to the measure fi^) if there exists an 
eigenvalue A e C so that for all g e Funct^ w (Vn), 

(5) E iQN.f}(q)W) = A ]T /(gjift). 

Using (4), we can rewrite this condition as (/, g) Dir = NX(f,g) e2 , which looks a 
good deal like the defining relation for cigenfunctions of A. In §3 we show how 
the eigenfunctions and eigenvalues of Qn with respect to the measure ^n relate 
to the usual notions of eigenvalues and eigenvectors. We show that all of the 
eigenvalues of Qn with respect to the measure fiN are positive in Proposition 4, and 
in Corollary 2 we prove that the eigenfunctions of Qn form a basis for Funct MN (Vn) ■ 
Let < Xi^n < ^2,n < ^3,n < ■ ■ ■ denote the eigenvalues of Qn, and write g^jv 
for the dimension of the eigenspace corresponding to the eigenvalue A^at. 

For each fixed i > 1 and TV > #X, let H N {i) C Funct^ (V N ) be an £ 2 - 
orthonormal basis of the eigenspace of Qn corresponding to the eigenvalue A^at- If 
the eigenspace is empty, set H.N{i) = (e.g., if TV < i). Define Tt(i) = Uat^aK*)- 
We think of the family H(i) as the set of all eigenfunctions of Qn corresponding to 
an ith eigenvalue. The family H(i) is obviously not unique, but we fix a choice of 
Ti.(i) for the remainder of the paper. 

Theorem 1 (Main Theorem). Fix i > 1. With the hypotheses and conventions as 
above, we have the following conclusions: 

(A) limAr^oo NXi y N = Ai(r). 

(B) There exists Nq — No (i) so that for all N > A^o, the dimension di,N of 
the eigenspace for Qn with respect to the measure [in corresponding to the 
eigenvalue A^at satisfies di t N — di(T). 

(C) The family H(i) is normal. The subsequential limits of H(i) lie in Zh Al (r) 
and contain an L 2 -orthonormal basis for the eigenspace of A corresponding 
to the eigenvalue Aj(r). 

As remarked above, we can even sharpen the statement of assertion (C) if the 
continuous part of the measure /i satisfies more smoothness properties. To reiter- 
ate, if [i — u{x)dx + J2 c j°~Pj an d w is C m away from the vertex set X, then the 
subsequential limits of the family H(i) are C m+2 on T \ X. 

The rate of convergence of eigenvalues and eigenfunctions is not studied in this 
paper. Our methods are purely existential, which is unfortunate due to certain 
interesting empirical data obtained by the students in the 2003 University of Georgia 
REU entitled "Analysis on Metrized Graphs" . The following two conjectures were 
made under the hypotheses that the models for T can be chosen with all edges 
having equal weights and that /i = dx: 

• There is a positive constant M such that |Ai(r) - iVAi,jv| < M ■ N~ 3 / 2 for 
all N. (The example in [KoFul] might lead one to believe that the error is 
more like 0(N^ 2 ).) 

• The scaled discrete eigenvalues AA^jv increase monotonically to the limit. 
This does not appear to be true for more general measures \i. 
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For more information about the REU and other data and conjectures resulting from 
it, see the official site: http://www.math.uga.edu/~mbaker/REU/REU.html 

In the next section we introduce an integral operator whose spectrum is in- 
timately related to the spectrum of Qn- We use it to show that Funct /JW (Vjv) 
admits a basis of eigenfunctions of the discrete Laplacian. In §4 we exhibit a reduc- 
tion of the Main Theorem which is technically simpler to prove. The proof of the 
Main Theorem will then proceed by induction on the eigenvalue index i (the case 
i = l will be identical to all others). We carry out the induction step in sections 
5-7. 

3. Integral operators and the spectral theory of Laplacians 

To begin, we recall the definition of the function j z (x, y). It is given by the unique 
continuous solution of A x j z (x,y) = S y (x) — 5 z (x) subject to the initial condition 
ic(C; y) — f° r au S/ST. Here A x denotes the action of the Laplacian with respect 
to the variable x. As the Laplacian of the j-function has no continuous part, one 
can show that j z {x, y) must be piecewise affine in x for fixed values of y, z. It is also 
true that j z (x,y) is symmetric in x and y, non- negative, jointly continuous in all 
three variable simultaneously, and uniformly bounded by 1. For elementary proofs 
of these facts see section 6 of [BF] . 

For the rest of this section, assume N > #A is a fixed integer. We define 
: T 2 -> C to be 

j^ N ( x ,y) = / j z {x,y)d/j, N (z) = 3q{x,y)m{<i)- 

Observe that j MJV (x, y) is also a piecewise affine function in x for fixed y. We require 
the following proposition, whose proof is given in [CR] as Lemma 2.16: 

Proposition 1. Let v be a signed Borel measure on T of finite total variation. 
There is a constant C v such that for each y G T, 

j v {x,y)dv{x) = C v . 

In light of this proposition, we define the kernel function g u : T 2 — > C to be 

g v {x,y) = j v {x,y) - C v . 

It follows that J r g v {x,y)dv{y) = 0. Also, the properties of the j-function men- 
tioned above imply that g v is symmetric in its two arguments and continuous on 
r 2 . Note that since T is compact, this forces g v to be uniformly continuous on T 2 . 

Now we recall the fundamental integral transform which effectively inverts the 
Laplacian on a metrized graph. For / € L 2 (T), define 

= j^9^{x,y)f{y)dy. 

It is proved in [BR] that ip^ is a compact Hermitian operator on L 2 (T) with the 
property that any element in the image of ip^ is orthogonal to the measure /i. In 
fact, we have the following important equivalence: 
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Theorem 2 ([BR], Theorem 12.1). A nonzero function f G Zh^(r) is an eigen- 
function of p^ with eigenvalue a > if and only if f is an eigenfunction of A with 
respect to the measure \i with eigenvalue A = 1/a > 0. 

By Proposition 15.1 of [BR] we know that any eigenfunction of A with respect 
to the measure fi lies in Zh M (r), and the above theorem allows us to conclude that 
eigenfunctions of have the same property. 

We now introduce a discrete version of the integral operator whose eigenvalues 
are related to the eigenvalues of the Laplacian matrix Qn in much the same way 
as in the above theorem. Define the discrete integral operator ipN ■ Funct(Viv) — > 
Funct(Viv) by 

<p N (h)(x) = J g^ N (x,y)h(y)dy N , where x G V N . 

The defining equation for p>N{h) only gives its values at the vertices, but recall that 
any function in Funct(Viv) is determined by its values on Vn by linear interpolation. 

The next lemma shows that the Laplacian matrix is essentially a left inverse for 
Pn, up to a scaling factor and a correction term. 

Proposition 2. If f G Funct(V}v), then for any q G Vn, we have 
{QnVn (/)} (q) = Jjf(<l) - f( x ) dx N^j m(q)- 

Proof. This is a restatement of Proposition 7.1 of [BR] (setting v = f(x)dxN, 
/i = f-i n 7 and using equation (4) to relate the Laplacian A to the discrete Laplacian 
matrix Qn)- □ 

Now we exhibit two useful facts about the kernel and image of the operator p>N- 

Lemma 1. Ker(<^jv) f] Funct^ (Vn) = 

Proof. Suppose / G Ker(t^Ar) P)Punct jUN (Vjv)- By Proposition 2, we find for each 
q G Vn that 

= {Qn^N (/)} (q) - ^f(q) - (J f{x) dx^j fi N {q). 

Put Cf = J T f(x)dxN- If Cf = 0, then / = 0. If Cf ^ 0, then the above equality 
implies f(q) = NCffj,N(q). But we then obtain the following contradiction to the 
fact that / G Funct^ w (Viv): 



L 



f(x)dn N (x) = NC f ]T M*) 2 ^ o. 

qieV N 

□ 



We remark that the previous proof actually shows the kernel of <pn acting on 
the space Funct(V^) consists of all scalar multiples of the piecewise afhne function 
defined by q i— > (J-N(q),q G Vat. 

Lemma 2. The operator pn is P-Hermitian with image in Funct^^V/v)- That 
is, if f G Funct(V/v), then 

J <p N (/) (x)dfi N {x) = 0. 



10 



X. W. C. Faber 



Proof. Suppose /, g G Funct(Vjv). As g flN is real and symmetric, we see 
(tPN(f),g) t 2 = (^J^g f , N (x,y)f(y)dy N ^j g{x)dx N 

= f(y) g^N (x, y)g{x)dx N ^j dy N 
= {f,VN{g)) P ■ 

Now recall that for any fixed y G T, the function g^ N (x,y) is orthogonal to the 
measure /xjv- Thus 

J^<p N (f){x)dn N (x) = ^g» N (x,y)f(y)dy N ^J dfi N (x) 

= J^f(y) (^J^g» N (x,y)diJ.N{x) s j dy N = o. 



□ 



It is interesting to pause for a moment to see how the eigenvalues and eigen- 
functions of Qn with respect to the measure hn relate to the usual notions of 
eigenvalues and eigenvectors. 

Proposition 3. The function f G Funct MAr (Vjv) is an eigenfunction for Qn with 
respect to the measure Hn if and only if there exists a constant A G C such that for 
all vertices q G Vjv , 

(6) {Q N f}(q) = A {/(g) - N [J f(x)dx N "j fi N {q) J . 

In particular, f is an eigenfunction for Q n with respect to the measure dx n if and 
only if f is a non-constant eigenfunction for Qn in the usual sense of a linear 
operator. 

Proof. The final claim follows from (6) because the integral term vanishes from 
the result when / is orthogonal to the measure dxN- 

If / G Funct MN (Vjv) is an eigenfunction for Qn with respect to the measure /Ujv, 
then there is A G C such that for all g G Funct MJV (Vat), 

E {QN.f}(q)W) = ^ E f(9)W)- 

We can rewrite this relation as TV (<3jv/ — A/, g) t2 = 0. Setting F — Qat/ — A/, we 
conclude that F is in the ^-orthogonal complement of the space Funct^ (Vn) (as 
a subspace of Funct(Viv))- One easily sees that the function q i— » /in(q) for q G Vn 
is a basis of the orthogonal complement. Thus F(q) = M/iN(q) for some constant 
M and all q G V N ■ 

The value {QNf}(q) can be interpreted as the weight of the point mass of A(/) 
at the point q using equation (4). As the Laplacian is always a measure of total 
mass zero, and the measure /in has total mass 1, we may sum the equation F(q) = 
M^ n {q) over all vertices q to see that 



M = -A E /(<?) = -^A / 



f(x)dx N . 



qev N Jr 
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This finishes the proof in one direction. The other direction is an immediate com- 
putation. □ 

Proposition 4. The eigenvalues oJQn acting on Punct(Vjv) are nonnegative. The 
kernel of Qn is 1- dimensional with basis the constant function with value 1. The 
eigenvalues of Qn with respect to the measure /Ujv are all positive. 

Proof. Using equation (4) and the notation for the weights on the edges of the 
model Gn in (3), we can expand the Dirichlet norm as 

N 

EiWKft) 5 qi (x) 



I/jvIIdit = J Mx) dA(f N ) = JJ N (x) 




/jvfe) {Qjv/at} (ft 



IMft)| 2 Y w ik~fN(qi) Y -M&Mj 

Note that wu = since Gn has no loop edges. Let us rearrange this last sum 
to be over the edges of Gat instead of over its vertices. Observe that summing 
over all pairs of vertices as above is equivalent to summing twice over all edges of 
the graph. We count w ik \f N (qi)\ 2 for each end of an edge ftft;. We also count 
—Wij /at (qi) /at (qj) and —WijfN(qi)fN(qj) for each edge qiqj. Here we have implic- 
itly taken advantage of the symmetry of Qn- This yields 

Wik \\jN\qi)\ + \jN(qk)' 2 

edges qiq k 



ll/jv||Dir= Y w ik{\fN(qi)\ 2 + |/jv(?fe)| 
edges qiq k 

(7) - Y ™y (Mft)Mft) + /jv(ft)/Ar(ft)) 



edges q^qj 

= Y w v - fN(qj)\ 2 ■ 

edges qiqj 

Note that in these sums we count edge qiqj = qjqi only once. 

If / S Funct(Viv) is an ^ 2 -normalized eigenfunction for Qn with respect to the 
measure fijy with eigenvalue A, then equation (7) and the defining relation for an 
eigenfunction shows that < H/Hdi,. = NX. Thus, A is nonnegative. If A = 0, then 
a more careful look at (7) shows that our eigenfunction / satisfies f(qi) = f(qj) 
for every pair of adjacent vertices qi,qj. Since Gjv is connected, we see that / is 
constant on the vertex set Vat. This proves the second assertion. 

Suppose further that / has constant value M. Then J r fdfiN = M. If / 6 
Funct^ N (Vat), we are forced to conclude that M = 0. That is, our function / is 
identically zero on Vn ■ This is horribly contrary to our definition of an eigenfunc- 
tion, and so we conclude that A > in this case. □ 

A function / e Funct(Viv) is called an eigenfunction for the operator ipx if there 
exists an eigenvalue a £ C such that fN(f) — ctf. Lemma 2 implies that any 
eigenfunction with nonzero eigenvalue must lie in Funct MJV (Vjv). Now we prove the 
theorem which relates the nonzero eigenvalues and eigenfunctions of the integral 
operator (p^ to the discrete Laplacian Qn- 
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Theorem 3. A nonzero function f £ Funct MJV (Vat) is an eigenfunction o/Qn (with 
respect to the measure ^n) with eigenvalue A if and only if it is an eigenfunction 
of with eigenvalue -^j. 

Proof. Suppose / is an eigenfunction of ipN in Funct MJV (Vat) with eigenvalue a. 
Proposition 2 implies that for each q £ Vn 

a{Q N f}(q) = {QNVN(f)}(q) = jjf(q) - ^f(x)dx N ^j p N (q). 

For g £ Funct MJV (Vjv), we may multiply this last equality by g(q) and sum over all 
vertices q £ Vn to see that 

a ]T {Q N f}(q)W) = jf E filW- 

q£V N q&V N 

The last sum vanishes because g is orthogonal the the measure /in- If a = 0, then 
/ = by Lemma 1; but eigenfunctions are not identically zero. Thus / satisfies the 
defining equation (5) for an eigenfunction of Qn with eigenvalue . 

Conversely, let / £ Funct MJV (Vjv) be an eigenfunction of Qn with eigenvalue A. 
Multiplying the result of Proposition 2 by 7VA and subtracting from the result in 
Proposition 3 shows 

QN{f-N\<p N (f)} = 0. 
By Proposition 4, the kernel of Qn is precisely the constant functions on T. Hence 
/ — N\(fiN(f) = M for some constant M. As / and <pn(/) are orthogonal to the 
measure /j,n, we know that integrating the equation / — N\(pn(J) — M against /kjv 
produces M = 0. Since all eigenvalues of Qn with respect to the measure ^n are 
nonzero, we deduce that ¥?jv(/) = j^t/- D 

Corollary 1. All of the eigenvalues of tfN acting on Funct^ (Vn) are positive. 

Proof. Apply the previous theorem and Proposition 4. □ 

Corollary 2. There exists a basis of Funct Mw (Vat) consisting of eigenfunctions of 
Qn with respect to the measure [in- 

Proof. The space Funct ([1 j V (Vjv) admits a basis of eigenfunctions of (fN by Lem- 
mas 2 and 1 and the finite dimensional spectral theorem. Now apply the previous 
theorem. □ 

4. Reduction to a weaker form of the main theorem 

Before embarking on the proof of the Main Theorem, we show the following 
reduction: 

Lemma 3 (Reduction Lemma). Suppose that for each fixed i > 1 the following 
assertions are true: 

(A) limjv^oo NXi^N = A;(r). 

(B') There exists No — No (i) so that for all N > No, the dimension of the 
eigenspace for Qn corresponding to the eigenvalue A^jv satisfies di.N < di(T). 

(C) The family H(i) is normal. The subsequential limits of H(i) lie in Zh /i (r), 
have unit L 2 -norm, and are eigenfunctions of IS. corresponding to the eigen- 
value A,(r). 
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Then Theorem 1 is true. 

The proof of the Reduction Lemma will require a few preliminary results, which 
will take up the majority of this section. 

Lemma 4. Suppose f G Zh(T) with the property that f is C 2 onT\X . Define fpj 
to be the unique function in Funct(Vjv) which agrees with f on the vertex set Vn ■ 
Then 

lim 1 1 ,/iV 1 1 Dir = H/llDir- 

Proof. We begin by recalling equation (7): 

(8) ll/iv||L= E ^l/ivfe)-/ivfe)| 2 . 

edges qiqj 

As T can be decomposed into a finite collection of segments {ejusing the preferred 
vertex set X, we can sum over the segments of T to get 

(9) ll/iv||L= E E m,\fN(q i )-fN(q 1 )\ 2 . 

segments e edges qtQj 

of r of Gn with Qi, qj G e 

The integral is linear, and there are only finitely many segments of T over which 
we wish to integrate, so it suffices to show that the inner sum in (9) converges to 
J e \ f'( x )\ 2 dx for any segment e of T (by equation (2)). Hence, we may assume for 
the remainder of this proof that T consists of exactly one segment e; that is, T is 
isometric to a closed interval of length 1 . 

The single segment e of T admits an isometric parametrization s e : [0, 1] — > e. 
The vertex set Vn corresponds to a partition = t\ < ti < ■ ■ ■ < tN = 1- We write 
f e = f o s e for case of notation. Now (8) takes the form 

JV-l 
N-l 

= e (*<+! - *i) 
i=i 

For each i there is t* e so that / e (t*+i) - / e (*i) = /e(*i)(**+i _ *») b y 

the mean value theorem. Hence 

AT-l 

11/^11^=2 **)l/e(*i)| 2 - 
i=l 

But this last expression is just a Riemann approximation to L \f' e (x)\ 2 dx. By 
Proposition 5.2 in [BR], we find that f' e is a continuous function on [0, 1] so that as 
N tends to infinity these approximations actually do limit to the desired integral. 

□ 

The following immediate corollary won't be of any use to us in our present task, 
but it is nice to include for completeness. 

Corollary 3. If f,g e Zh(T) 7 then we have the limit 

lim (/iv,5Jv) Dir — > (/,5> Dir , 



fe(U+l) - fejtj) 

ti+i — u 



14 



X. W. C. Faber 



where /n and gN are affine approximations of f and g as in the statement of 
Lemma 4. 

Proof. This is an easy consequence of Lemma 4 and the polarization identity 

3 



(fN,g N ) Dir = \^i n \\fN + i n g N \\l h 



4 

which may be checked easily by expanding the norms on the right-hand side. Here 
i denotes a fixed complex root of —1. □ 

Here is a technical lemma that will be needed in the Approximation Lemma 
(Lemma 6) and again in later sections. 

Lemma 5. Suppose for each N that is a function in Punct(Vjv) such that 
the sequence {kjy} converges uniformly to a (continuous) function k : T — > C. 
Then ||/W||i 2 — > II^IIl 2 as N — > oo. If {k' N } is another such sequence of functions 
converging uniformly to some k' : T — > C, then (fcjv, k' N )p — > (k, k') L2 as N — > oo. 

Proof. Suppose N is so large that fe^ is uniformly within e of k. On one hand, we 
have 

|2 _ f \u__<„\\2. 



I^nWp = I \k N (x)\ dx N 



L 

j {\k(x)\ + e) 2 dx N 



< 

'i' 

|| 2 2 +2 e ||fc||, 1+£ 2 

2 i o-III.ll _ i ~2 



^||fc||£ 2 +2 £ ||fc|| il+ ^ 

The convergence in the final step follows from weak convergence of dxM to dx. 
Letting e — > shows that limsup^^^ ||fcjv||< 2 < II^IIl 2 - 



Similarly, 



\kpf(x)\' 2 dxN > / (\k{x)\ ~ e) 2 dx n 



| 2 2 -2e||fc||, 1 +e 2 

2 o-ll7.ll . i -2 



— ►||fc||£ 2 -2 e ||fc|| il + e '. 

Now we see that ||fc||_L 2 < liminfjv_»oo ||fcjv||^ 2 - 

The final statement follows by the polarization identity as in the proof of Corol- 
lary 3. □ 

Lemma 6 (Approximation Lemma). Fix m > 1 and suppose f £ Zh M (r) is an 
L 2 -normalized eigenfunction of A with corresponding eigenvalue X m (T). For each 
N, let /at be the unique function in Funct(Vjv) which agrees with f at the vertices 
in Vjv. Assume also that assertions (A), (B'), and (C) of the Reduction Lemma 
(Lemma 3) are satisfied for j < m — 1. Given any subsequence of {Gn}, there exists 
a further subsequence such that for each e > there is a positive integer N\ with 
the property that for every N > N\ with Gn in our sub-subsequence, we can find a 
function /at € Funct MJV (Vn) for which the following conditions hold simultaneously: 

(i) WjNy = 1; 

(ii) /at is £ 2 -orthogonal to every h G HnU), j = 1, • • • , m — 1; and 
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(iii) H/ivllDir - 11/jvllDir 



< e. 



Proof. As condition (B') of the Reduction Lemma holds, and each dj(T) is finite, 
we may pass to a subsequence of models {GV} such that dj^ is independent of N 
for every j < m — 1 and all TV sufficiently large in the subsequence. This implies 
that for large N, the set Tv = UjLi 1 ^n{j) is finite with cardinality independent of 
N. Denote this cardinality by T, and let us label the elements of Tv as h} N , . . . , hjj. 
The set T/v is an £ 2 -orthonormal system by definition of the sets HnU) an d the 
fact that eigenfunctions corresponding to distinct eigenvalues are £ 2 -orthogonal (a 
consequence of the definition of cigenfunction) . 

Write Bm = J r fd/iN', we will abuse notation in what follows and let Bm also 
denote the constant function on T with value Bjy. Define 



In - Ej=i (/jv, h 3 N ) h J N - B N 



N 



In - Ej=i {In, h J N j p h 3 N - B N 



It will be shown in a moment that for TV sufficiently large lying in a certain subse- 
quence, the denominator of this expression is nonzero. 

The definition of H(j) implies that each h J N is orthogonal to the measure /xjv, so 
a small calculation shows Jn lies in Funct^ N (Vn)- Evidently condition (i) holds. 
Also, we note that Bn is an eigenfunction of ip^ acting on Funct(Vjv) with cor- 
responding eigenvalue zero. Any two eigenfunctions of a self-adjoint operator that 
have distinct associated eigenvalues are orthogonal, and hence Bm must be £ 2 - 
orthogonal to each h J N . Another quick calculation shows that h 3 N is orthogonal to 
/jv for all j < m — 1, which is condition (ii). 

It remains to check that condition (iii) holds for our choice of /at upon passage to 
a further subsequence. First observe that the definition of an eigenfunction for Qjy 
implies that distinct elements of T/v are orthogonal with respect to the Dirichlet 
inner product. Also, (S7v,5) Dir = for all g € Funct(Vjv) since constant functions 
arc in the kernel of Qm- Expanding the Dirichlet norm and simplifying using these 
two observations gives 



T 



HL/VllL-E^w {fN,h% 

3=1 



T 

2 



(2 



Dir 

where Q.n is the eigenvalue associated to the eigenfunction h J N . Of course, Q.n — 
\k,N for some k < m — 1. 

We have T sequences of functions {h]^}, . . . , {hjj}, each of which lies in a par- 
ticular H(j). The normality of each H(j) allows us to find a further subsequence 
of models and a set of functions T = {h 1 , . . . , h T } so that h 3 N — > h? uniformly on 
T along this subsequence. For the rest of the proof we will assume that N lies in 
the subsequence we have just specified. 

As T/v forms an ^ 2 -orthonormal set for each N, we see that T forms an L 2 - 
orthonormal set of eigenfunctions for A (assertion (C)). The associated eigenvalues 
of these functions are all strictly smaller than A m (r), and hence each h 3 is T 2 - 
orthogonal to our initial function /. Lemma 5 implies that the inner product 
(^fN-,h 3 N ^j ^ tends to (/, h J ) L2 = as N tends to infinity. 
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By Lemma 4, weak convergence of {dx^} to dx, the previous paragraph, and the 
hypothesis that the scaled eigenvalues NXj^n converge to Aj(T) for each j < to — 1, 
we conclude that 



(io) II/aHIL-E^." (-^ 

as iV tends to infinity. 



i 2 



lr>ir 



E<v </•/'''> 



L 2 I 



llJir 



We also note that the function X) 7 =i \fN,h 3 N ) h 3 N + Bn tends to zero uni- 



i 2 



formly. Here Bn — > by weak convergence of {/i at} to \i. The Minkowski inequality 
and Lemma 5 imply that 



fa - E \f^,K 



7 ^AT 



< II/a/II^ + 



L 2 +o(l), 



where the error term depends only on N. The Minkowski inequality also gives an 
identical lower bound so that 



(11) 



3 = 1 



N 



= 1 



as TV goes to infinity. This shows that the denominator of fa is in fact nonzero for 
TV sufficiently large. 

Finally, we use the convergence relations (10) and (11) along with Lemma 4 to 
conclude that as TV tends to infinity 

„ H/iv|lD ir -Ej=1^0,iv|(/iV,^ 12 



N 



II/. 



t 2 



/Ar-E;=i(/iv,^), h%-B 



II/. 



ATllDi) 



0. 



This completes the proof of condition (hi) 



□ 



And now we require one final proposition which gives us a characterization of 
the rath largest eigenvalue of Qn as a minimum of the Dirichlet norm on the £ 2 
unit circle. 

Lemma 7. Fix N > #X. For a given to > 1, let fa — \Sj=i 7~Ln{j)- Denote by 
the £ 2 -orthogonal complement of the span of fa inside the space Funct MAr (V/v). 
Provided that ^ 0, we have 



NX m ,N = min 



\9\\h 



gefa \\g\\% 



Proof. As Funct AIN (Fjv) admits a basis of eigenfunctions of Qn (Corollary 2), 
there exists an £ 2 -orthonormal basis {/i, • • • , / r } for . Let us assume that the 
corresponding eigenvalues are 71 < 72 < ■ ■ ■ < 7r- Note that 71 = \ m ,N since 
consists of all of the eigenfunctions of Qn with eigenvalue strictly larger than 
A m -i,jv- 
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To prove the proposition it suffices to consider the minimum over all g G Tjq with 
unit £ 2 -norm. Using our basis, write g = Y^j=i a jfji where the complex numbers 
dj satisfy ^\ \dj\ 2 = 1. Since (fj,fj) uil = N-fj, we see that 

r r r 

M Dir - E n 2 /i>Dir = Ar Ei ffl i^ ^ ^ E m 2 = NX ^- 

Note that if we take g G Funct AlAr (Viv) to be an £ 2 -normalized eigenfunction of Qn 
with eigenvalue X m .N, then equality is actually achieved in this last computation. 

□ 

Finally, we may now return to the 

Proof of Reduction Lemma. We assume that conditions (A), (B'), and (C) 
hold for every i > 1. Clearly we only need to check that assertions (B) and (C) in 
the Main Theorem are true. If assertion (B) fails for some value of i, then we may 
choose a minimum such i. There exists a positive integer D strictly smaller than 
di(T) and an infinite subsequence of models {Gn} with d^N = D for all Gn m 
the subsequence. For each such N, select distinct functions h} N , . . . , G TLnH)- 
Using the normality of H(i) as given by hypothesis (C), we may pass to a further 
subsequence and assume that there are D limit functions h 1 , . . . ,h D , with h 3 N — > /i J 
uniformly along this subsequence for each j. Again using assertion (C) we find that 
each hP is an eigenfunction of A with associated eigenvalue Aj(r). As D < di(T), 
there is some eigenfunction of A with eigenvalue Aj(r) which is not in the span of 
{ft 1 , . . . , h D }. Choose such a function with unit L 2 -norm and denote it by /. 

Now choose a sub-subsequence for which the conclusion of the Approximation 
Lemma is satisfied for the function / in the case m = i + 1. Given e > and N 
large, we pick /jv as in the Approximation Lemma. We apply Lemma 7 in the case 
m = i + 1 to see that 

N\i + i n = min < n/^112 < \\f N \\l ii + s. 

We remark that ^ for N sufficiently large since its dimension is at least TV — 
1 — Y^ij=i by assertion (B'). Letting N go to infinity through our subsequence 

and using assertion (A) and Lemma 4, we find that 

wn < 11/112*+ e. 

We know that (/, ff) Dir = Aj(r) (/, g) L2 for all g G Zh AI (r) since / is an eigenfunction 
of A. In particular, setting g = f tells us that ||/||D ir = Aj(r). Taking e sufficiently 
small provides a contradiction since Aj(r) < Aj + i(r). Thus assertion (B) must be 
true for all i. 

Now we prove assertion (C). As assertion (B) holds for each fixed i, we know 
that Hjv(i) consists of d = di(T) pairwise ^-orthogonal functions on T when N is 
large. Let h} N , . . . , h% be the elements in Wjv(i). By normality of H(i), we can pass 
to a subsequence such that each of these d sequences {ft^} converges. Now use 
Lemma 5 to see that the limit functions must be pairwise L 2 -orthogonal and have 
unit L 2 -norm. The dimension of the eigenspace of A corresponding to Aj(r) is d, 
so the subsequential limits just constructed form a basis for this eigenspace. □ 
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5. Normality of the family 7i(i) 

For the remainder of the paper we assume that assertions (A), (B'), and (C) of 
the Reduction Lemma hold up to the case i — 1, and we prove that they are true 
for the case i. 

Theorem 4 (Eigenvalue Convergence: Part 1). We have 

lim sup NX^n < Aj(T). 

Proof. Suppose the result is false and pick a subsequence of models and a positive 
constant 5 so that NXi t N > Xi(T) + 5 for all N sufficiently large in our subsequence. 
Let / be an L 2 -normalized eigenfunction for A in Zh^(r) with eigenvalue Aj(r). 
We now select a further subsequence as in the Approximation Lemma of §4. For 
any e > and TV sufficiently large in our sub-subsequence, we can find Jn and 
apply the Approximation Lemma and Lemma 7 to get 

Dir , || f || 2 , _ 

X,{i ) + d < JMX ifN = mm < < ||/jv|| D ir + e - 

geT^ \\9\\p \\fN\\ 2 e 2 

We remark that T^- is nonempty for N sufficiently large by assertion (£?') (which we 
are assuming holds up to i — 1). Letting TV tend to infinity through our subsequence 
and applying Lemma 4, we arrive at the statement 

Ai(r) + <y < ll/HL + e, foranye>0. 

As H/llrjir = -^i(r), taking e sufficiently small produces a contradiction. □ 

Corollary 4. The family TL{i) has uniformly bounded Dirichlet norms. 

Proof. By definition of an eigenfunction of Qn with respect to the measure /xjv, we 
have immediately that ||/i||D ir = NXi^jy for any h E Wjv(i). The previous theorem 
shows that ||/i||j)i r cannot be arbitrarily large for h lying in H(i). □ 

Theorem 5. (a) The family H(i) is equicontinuous on V. 

(b) TL{i) is uniformly bounded on T. 

(c) H(i) is a normal family. 

Proof, (a) Suppose M is a positive real number such that ||/i||Dk < M for all 
h E TL{i). Now suppose that e > is given and that we have any pair x, y E F 
with dist(x,y) < e 2 /M 2 , where dist(-, •) is the metric on T. Let 7 C T be a 
unit speed path from x to y. The fundamental theorem of calculus and the 
Cauchy-Schwarz inequality imply that 



\h{x)-h{y)\ 



f h'(t)dt 

< (y <ttj 1 (j \h\t)fdt 



< dist(z,y) 1/2 ||/i||Dir < £• 

The integrals in the above estimate should be viewed as path integrals in which 
we have implicitly chosen compatible orientations for consecutive segments 
along the path 7. 
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(b) It suffices to show that the family H(i) is uniformly bounded on [j N Vn as 
this constitutes a dense subset of V and all of our functions are continuous. 



The last inequality follows because the distance between any two points in 
r is at most 1. Now Corollary 4 and the assumption that the measures /in 
have uniformly bounded total variation show that there is a constant M, 
independent of N, so that \h(v)\ < M for any choice of h E Hn and v E Vn- 
Take any h E fi,(i) now. Choose N so that h E Wjv(*)- Note that if x G T, 
then x lies on some segment e of Y with respect to the vertex set Vn- By our 
choice of N, we find h is affine on e. Taking y,z EVn to be the endpoints of e, 
we know that there is t E [0,1] so that h(x) = (l—t)h(y) + th(z). The triangle 
inequality implies \h(x)\ < max{|/i(t/)|, \h(z)\} < max Vie v N ^ M. 

(c) Apply the Arzela-Ascoli theorem using parts (a) and (b) to satisfy the neces- 
sary hypotheses. 



6. Convergence and approximation results 

This section contains the proofs of various technical lemmas regarding the inte- 
gral operator ip^ and its relation to the discrete integral operator ip n- 

Lemma 8 ([Rud], p. 168, Exercise 16). Suppose {k^} is an equicontinuous sequence 
of functions on the metrized graph T such that {k^} converges pointwise. Then 
{k N } converges uniformly on T. 

Proof. Fix e > 0. By equicontinuity there exists 5 > such that \kN(x) — fcjv(y)| < 
e/3 for every N and x, y G T with dist(x, y) < S. By compactness, we can cover V 
with a finite number of balls of radius 5/2. Select y\, . . . ,y m e T so that at least 
one yi lies in each ball. 

For j = 1, . . . , m, let Aj be a positive real number such that M,N > Aj implies 
\^N(Uj) — kM{yj)\ < e/3. This is possible because {fcjv} is pointwise convergent. 
Set A = max{ A,- : j — 1, . . . , to}. 

Suppose M, N > A and that x E T. Pick yt so that dist(x, yi) < 5. Then 

\k N (x) - k M (x)\ < \k N (x) - k N (yi)\ + \k N (yi) - k M {yi)\ + \k M (Vi) - k M (x)\ < e. 

Note that our choice of M, N did not depend on x. This shows {k^ is uniformly 
Cauchy. □ 

Lemma 9. Suppose that {kj\[} is a sequence of continuous functions on T which 
converges uniformly to a function k : T — > C. Then (^(/cn) ~ > uniformly 
with N. 



Recall that if h E Hat(i), then J r hdfijy = 0. By the proof of part (a) it 
follows that for any vertex v E Vn and h E Wat(«), 




< ||A|| D irM(r). 



□ 



20 



X. W. C. Faber 



Proof. Our strategy will be to use Lemma 8. Fix x e Y. Then, since g^(x,y) 
is continuous on T 2 , it must be uniformly bounded (compactness). Also, since 
{/sat} converges uniformly to a continuous function k, we find that fcjv(y) is uni- 
formly bounded for all N, y. Thus Lebesgue dominated convergence guarantees 
that <^(fc/v) — ► <fin(k) pointwise. 

Suppose that M is a positive real number such that |fcjv(j/)| < M for all N,y. 
Also, since g^ is continuous on a compact space, it must be uniformly continuous. 
Thus, given e > 0, we can select S > such that \g^(x,y) — g^,(x',y)\ < e/M 
whenever x,x',y G T with dist(x,x') < 5. Hence for all x,x',y <G T such that 
dist(x,x') < S, we have 

\fn{kN)(x) - <Pp(k N )(x')\ < J \gp(x,y)-g li (x',y)\-\k N (y)\dy <e. 

This shows the desired equicontinuity. □ 

Lemma 10. As N — ► oo 7 we have g^ N — > g M uniformly on T 2 . 

Proof. We intend to apply Lemma 8. 

First we show pointwise convergence. By definition of g^ N and g^, it suffices to 
show 

• Jun ( x > V) -> 3n( x i V) f° r eacn P air x, y <G T, and 

• J T 3w(x,y)diJ, N (x) -> J r jn(x,y)d/j l (x) for all i/GT. 

The former assertion follows by weak convergence of measures. As for the second, 
fix e > and a basepoint yo € T. As j z (x,y) is a continuous function on T 3 (a 
compact space), it is uniformly continuous. So there exists <5 > such that for any 
x,x',z £ T with dist(x,x') < S, we have \j z (x,yo) — j z (x',yo)\ < e. Now for such 
x, x', it follows that 

l-WO^yo) - jn N (x',yo)\ < J \jz{x,yo) - jz(x' ,y )\d\n N \(z) <e\n N \(T). 

As the total variations |/njv|(r) are uniformly bounded independent of N, we see 
{Jun ( x , yo)} is equicontinuous in the variable x. By Lemma 8 we find j^ N (x, yo) — > 
jn(x,yo) uniformly in x. For ease of notation, set I^(ct) = J r j u (x,y )da(x) for any 
pair of measures v, er. Now suppose that N is so large that j^ N {x, yo) is uniformly 
within e of j^(x,yo), and that I^Hn) is within e of I^(^i) (by weak convergence). 
Then 

\I» n (Vn) - ^(/i)| < \I» N (fJ-N) - ^(mjv)| + I^(Mjv) - 
(12) < J^\j» N (x,yo)-j»(x,y )\ d\fi N \(x) + |/ m (a6jv) - 

<e| M jv|(r)+e. 

By Proposition 1, we see that J r j llN {x,y)dnN{x) and J r j ll (x,y)d/i are indepen- 
dent of As e was arbitrary, we conclude from (12) that J r j^ N (x,y)dfiN{x) — > 
Jr3n( x , y)dfi pointwise for all y G T. 

All of this allows us to conclude that g^ N {x,y) — ► g fl (x,y) pointwise on T 2 . It 
remains to prove {y MN (x, y)} is equicontinuous jointly in x and y. By uniform 
continuity of j on T 3 , we know there exists S > such that for all x, x', y, y', z with 
dist ((x, y), (x',y')) < 5 (some metric on the product space T 2 ), it must be that 
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\jz(x,y) -j z (x',y')\ < e. Now for any x,x',y,y' G T with dist ((x, y), (x' , y'j) < 5, 
we have 

\9m(x,y) - 9m(x',y')\ = \j»N{x,y)-jn N {x',y')\ 

< J \j z (x,y)-j z (x',y')\d\» N \(z) < s\f, N \(T). 

This proves the equicontinuity of {g^ N {x,y)}. □ 

Lemma 11. Given any convergent sequence {Hn} C H(i) with limit function h 
such that Km G Hat(«) for each N, we find that (fiN(hN) — ► <Pn{h) uniformly on T. 

Proof. We use Lemma 8 again. First we show pointwise convergence. Fix x E T 
and e > 0. As <? MJV — > g M uniformly and hjy ^ h uniformly, we can suppose 
that N is large enough to guarantee g^ N hfq is uniformly within e of g^h for all 
x, y G T. We can also suppose that N is so large that J r y)h(y)dyN is within 
e of J r g^(x,y)h(y)dy. Now we have 



j 9n N (x,y)h N (y)dy N - J g fl (x,y)h(y)dy 

9 UN 

(x, y)h N (y)dy N - J g^x, y)h(y)dy N 

(x, y)h{y)dy N - J g^x, y)h{y)dy 



< 



9n\ 

T 

< j^\9^ N {x,y)h N {y) - g^{x,y)h{y)\dy N 

< 2e. 



Hence pointwise convergence holds. 

Now we wish to show equicontinuity of the sequence of functions {</?iv(ftjv)}. The 
sequence {g^ N (x, y)} is equicontinuous on T 2 by the proof of Lemma 10, and the 
sequence {Hn} is uniformly bounded by some positive constant M. Thus we have 

\(p N {h N }(x) - f N {h N }(x')\ < M J \g^ N {x,y) - g t _ tN (x' ,y)\dy N , 

and the integrand can be made arbitrarily small by taking x close to a;'. □ 

Lemma 12. Suppose that {h^} C T~t(i) is any convergent sequence with limit 
function h such that ftjv G Wjv(i) for each N . For each e > there exists a real 
number N 2 such that for all N > N 2 , we have simultaneously 

• H'MMIk 2 > II Vjv (ft at) ||^ -e; 

• \\h\\ L 2 < \\h N \\ £ 2 + e. 

Proof. Observe that ft at — > ft uniformly, and so y>jv(ftjv) —> fn{h) uniformly by 
Lemma 11. Now use Lemma 5 to assert that for all N large, ||ftjv||f2 is within e of 
||ft||i2 and \\<PN(hN)\\e 2 is within e of \\<p^(h)\\ L 2 . □ 

Lemma 13. Suppose ft : T — > C is a continuous function such that J r hdu = 
and ft is L 2 -orthogonal to all of the eigenfunctions of A with associated eigenvalues 
Ai(r), . . . , Aj_i (r). For each e > 0, there exists ft G Zh M (r) that satisfies the 
following two conditions: 
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(i) h is L 2 -orthogonal to all of the eigenfunctions of A with associated eigenvalues 
Ai(r), . . . , Ai_i(r); and 

(■■\ \\Vv(h)\\ L 2 \hf»(h)\\ L 2-e 

y > ii r,i „ 



llL 2 „a\\ L 2+e 

Proof. Using the vertex set X , we may decompose T into a finite number of seg- 
ments {ei, . . . , e r }, each isometric to a closed interval. As h is continuous, we can 
apply the Stone- Weierstrass theorem to each segment of L to get uniform polyno- 
mial approximations of h\ ej with as great an accuracy as we desire. Moreover, we 
may define our approximations so that their values agree with the values of h at 
the endpoints of each segment, and hence we may glue our approximations together 
to get a uniform approximation of h. Evidently such approximations must lie in 

zh(r). 

Let us perform this procedure to construct a sequence of functions {/„} C Zh(r) 
that converges uniformly to h. Let £1, . . . , £ p be eigenfunctions of the Laplacian A 
that form an L 2 -orthonormal basis for the direct sum of the eigenspaces associated 
with the eigenvalues Ai(T), . . . , Aj-i(r). Now set 

p , 

H n {x) = f n (x) - (fn,^j) L 2 £,j(x) - / f n dfl. 

3 = 1 Jr 

As each ^ is orthogonal to constant functions (constants are eigenfunctions of ip^ 
with eigenvalue 0), we see that H n £ Zh M (r) for all n and condition (i) of the 
lemma is satisfied for all H n . Since /„ tends uniformly to h, we see that the inner 
products {f n ,£,j) L 2 tend to zero as n tends to infinity. By dominated convergence, 
the integral J r f n d\x tends to zero as well. Thus H n converges uniformly to h on L. 

By Lemma 9, we know that ip^,(H n ) — ► <f^(h) uniformly. Since the L 2 -norm 
respects this convergence, for any e > we may take n sufficiently large and set 
h = H n to obtain condition (ii). □ 

The eigenvalues of ip^ constitute a sequence of positive numbers tending to zero. 
If we label these distinct values by ai(r) > 0:2(1") > 0:3 (r) > • • • , then o^(L) = 
1/Aj(r) for all i (see Theorem 2). We have the following classical characterization 
of the ith eigenvalue of ip^: 

Proposition 5. Let S(i) denote the span of the eigenfunctions of A with respect 
to the measure /i that are associated to the eigenvalues Ai(r), . . . , Aj_i(r). If S(i) 1 - 
is the L 2 -orthogonal complement of this space inside Zh M (r), then we find that 

(13) ai{T) = sup — f— . 

Moreover, there exists an eigenfunction F £ S(i) which realizes this supremum. 

Proof. In section 8 of [BR], it is proved that ip^ is a compact Hermitian operator 
on L 2 (T); hence, by the spectral theorem, the eigenfunctions of cp^ form an or- 
thonormal basis for L 2 (T). In section 15 of [BR] it is shown that the eigenfunctions 
of must lie in Zh M (r). We may now identify <S(z)^ with a subspace W of L 2 (T) 
that is orthogonal to all of the eigenfunctions of A associated to the eigenvalues 
Ai(r), . . . , Aj_i(r); moreover, we see that ip^ maps this subspace into itself. 
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The space Zh(T) is dense in L 2 (F) because, for example, smooth functions are 
dense in L 2 . Thus the supremum on the right side of equation (13) is nothing 
more than the operator norm of (p^lvK- By [Yng] Theorem 8.10, we find that either 
IIVmIwHI or — HwdVII is an eigenvalue of ip^ acting on W . The nonzero eigenvalues 
of ifp are all positive, so we may eliminate the latter possibility. We conclude 
that there exists an eigenfunction F E W with associated eigenvalue a for which 
llwlwll — ^(^OII^/II-^IU 2 = a - Evidently a = cci(r) since the operator norm is 
defined as a supremum and oti (T) is the largest eigenvalue of ip^ acting on S(i) x . □ 



7. Proofs of assertions (A), (B'), and (C) 

Let cti.N > C(2,n > «3,Ar > • ■ ■ denote the nonzero eigenvalues of the operator 
<Pn described in §3. There we saw that for each i, the eigenvalues of (fN are related 
to the eigenvalues of Qm with respect to the measure /in via iVA^jv = l/a^jv. 

Recall that we are assuming assertions (A), (B'), and (C) of the Reduction 
Lemma hold up to the case i — 1. 

Theorem 6 (Eigenvalue Convergence: Part 2). We have 

Xi(T) < liminf NX hN . 

AT— >oo 

Proof. Suppose the theorem is false. As N\ iy jy — l/ a i.N and Aj(r) = l/a^T), our 
supposition is equivalent to aj(r) < limsupjy^^ a^jv- We may pick a subsequence 
of {Gn} and a positive constant r\ so that c^(r) + 77 < a itN for each N in the 
subsequence. For each such N, we also pick Hn G Wjv(i)- By normality of H(i), 
upon passage to a further subsequence we may assume that {/ijv} tends uniformly 
to some continuous function h as N goes to infinity through this subsequence. 

Fix e > and select h as in Lemma 13. For N large, the conclusions of Lemma 12 
are true. It now follows that for any N sufficiently huge in our subsequence, we 
have 

ai (T)= sup y ^ F , )h2 , by Lemma 5 

> \Vp£^v L shlce ~ heS(i) ± 

~ \\h\\v 

> " y ,ff llL r g . by Lemma 13 
\\h\\ L 2+e 

^ll^yilf- 26 , by Lemma 12 

since h n is an eigenfunction of ipN 



> 



\\h\\ L ^ 


- £ 


\\<PnQi n )\ 


\p - 2e 




+ 2e 


Oti,N /iJV 


£2 - 2e 


\\h N \\p 


+ 2e ' 


ati^N - 2£ 




l + 2e 






- 2e 



l + 2e 

This furnishes us with a contradiction when e is sufficiently small. □ 
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Clearly the amalgamation of Theorems 4 and 6 prove that assertion (A) of our 
Main Theorem holds for i. As for assertion (B') and the remainder of assertion 
(C) in the Reduction Lemma, we have the following two corollaries. 

Corollary 5. The sub sequential limits ofTL(i) are eigenf unctions of A with asso- 
ciated eigenvalue Ai(r). The limits have unit L 2 -norm and lie in Zh M (r). 

Proof. Let {h^} be a convergent subsequence of Tt(i) with limit function h. We 
know that (pN(hpf)(x) = - Hn{x) for all x e F. Now let N go to infinity and 

apply Lemma 11 and assertion (A) to see that ^p^{h) = x\r) ^" Hence h is an 
eigenfunction of ip^. By Theorem 12.1 and Proposition 15.1 in [BR] we see that h 
is an eigenfunction of A lying in Zh M (r). Lemma 5 shows ||/i||z,2 = 1. □ 

Corollary 6. For N sufficiently large, it is true that c^jv < di(T). 

Proof. Set d = di(T). Suppose that that di^ > d for some infinite subsequence of 
{Gn}- For each such N, select distinct functions h]^, . . . , h^ 1 e Hn^)- Using the 
normality of H(i), we may pass to a further subsequence and assume that there are 
d + 1 limit functions h 1 ,. .., h +1 , with h 3 N — > h? uniformly along this subsequence 
for each j. By Corollary 5, we find that each h? is an eigenfunction of A with 
associated eigenvalue Ai(T). But the functions h x N , . . . , hff 1 are € 2 -orthogonal for 
each N, and Lemma 5 implies that the limit functions must be L 2 -orthogonal. As 
the eigenspace corresponding to the eigenvalue Ai(T) has dimension d, this is an 
absurdity. □ 

We have now exhibited the induction step in our proof of the Main Theorem; 
i.e., the proof is complete. 

8. Acknowledgements 

This work could not have been completed without the gracious support of the 
NSF VIGRE grant at the University of Georgia or the support of the Columbia 
University mathematics department. I would like to thank Matthew Baker and 
Robert Rumely for introducing me to the subject and for their most helpful guidance 
toward completing this project. I would also like to thank David Glickenstein for 
pointing me toward the relevant work being done in China and Japan. The members 
of the 2003 UGA REU on metrized graphs deserve a round of applause for their hard 
work and enthusiasm. My wife has always provided me with a loving environment 
in which to work, and for this she should be praised most highly. 



References 

[BF] Matt Baker and Xander Faber. Metrized graphs, electrical networks, and Fourier analysis. 
Submitted. 

[Bo] Bela Bollobas, Modern Graph Theory, Springer- Verlag, New York, 1998. 

[BR] Matthew Baker and Robert Rumely. Harmonic analysis on metrized graphs. To appear in 
Canad. J. Math. 

[CR] Ted Chinburg and Robert Rumely. The capacity pairing. J. reine angew. Math. 434 (1993), 
1-44. 



Spectral Convergence 



25 



[KoFul] Koji Fujiwara. Convergence of the eigenvalues of Laplacians in a class of finite graphs. 

Geometry of the spectrum. (Seattle, WA, 1993) Contemp. Math., 173, Amer. Math. Soc, 

Providence, RI, 1994. 115-120. 
[KoFu2] Koji Fujiwara. Eigenvalues of Laplacians on a closed Riemannian manifold and its nets. 

Proc. Amer. Math. Soc. 123 (1995), no. 8, 2585-2594. 
[KcFu] Kenji Fukaya. Collapsing of Riemannian manifolds and eigenvalues of Laplace operator. 

Invent. Math. 87 (1987), no. 3, 517-547. 
[Ku] Peter Kuchment. Quantum graphs: I. Some basic structures. Waves Random Media. 14 

(2004), S107-S128. 

[Mo] Bojan Mohar. The Laplacian Spectrum of Graphs. Graph Theory, Combinatorics, and Ap- 
plications. Vol. 2, Wiley, 1991, 871-898. 

[Rud] Walter Rudin. Principles of Mathematical Analysis. 3rd cd., McGraw-Hill, 1976, New York. 

[Yng] Nicholas Young. An Introduction to Hilbert Space. Cambridge University Press, 1988, New 
York. 

Department of Mathematics, Columbia University, New York, NY 10027, USA 
xander@math.columbia.edu 



